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Abstract 

We investigate scalar and spinor field theories in a constant mag- 
netic field at finite temperature and chemical potential. In an external 
constant magnetic field the exact solution of the two-point Green func- 
tions are obtained by using the Fock-Schwinger proper-time formalism. 
We extend it to the thermal field theory and find the expressions of 
the Green functions exactly for the temperature, the chemical poten- 
tial and the magnetic field. For practical calculations the contour of 
the proper-time integral is carefully selected. The physical contour is 
discussed in a constant magnetic field at finite temperature and chem- 
ical potential. As an example, behavior of the vacuum self-energy is 
numerically evaluated for the free scalar and spinor fields. 
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1 Introduction 



Some of the interesting cosmological and astrophysical situations are found in 
the state with high density, temperature and strong magnetic field. Neutron 
star is a dense object which has a large chemical potential. Recently it is ob- 
served that some of the neutron stars have extremely strong magnetic field. ""^ 
A primordial magnetic field in the early universe is also interesting.^''^ To un- 
derstand the physics of such situation we consider the quantum field theory 
in an external magnetic field at finite chemical potential and temperature. 

One of the fundamental objects in the quantum field theory is a two-point 
Green function. The exact form of the Green function is necessary to deal 
with a strong magnetic field. Much interest has been paid to obtain the 
Green function under external fields. Schwinger found the exact expression 
of the Green function in an external magnetic field by using the proper-time 
formalism in 1951.^ The proper-time method is extended to deal with the 
thermal system in Ref. 5. Trace of the Green function corresponds to the 
vacuum self-energy of the free field. It is obtained by QED effective action 
and evaluated in a constant electromagnetic field at finite temperature^''' and 
at finite chemical potential.^'^ Variety of approaches are used to discuss the 
contributions from both the temperature and the chemical potential in Refs. 
10-13. 

In the present paper the proper-time formalism is re-considered in the 
imaginary time form of the thermal field theory. We modify the formalism to 
introduce both the temperature and the chemical potential exactly. In most 
of previous analysis the proper-time integral was analytically performed by 
the Landau level expansion. Since results are analytically continued to the 
wide range of parameters in the Landau level approach, these results may 
contain some approximation for the combined effect of temperature, chemical 
potential and external magnetic field. To obtain physical results we carefully 
choose the contour of the proper-time integral. Here the explicit form of 
the scalar and fermion Green functions is written down in a proper-time 
form to discuss the physical contour in a constant magnetic field H at finite 
temperature T and chemical potential /x. As an example, we numerically 
calculate the vacuum self-energy for the free scalar and fermion fields and 
discuss the contour dependence of the proper-time integral. 
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2 Scalar Two-point Function at Finite T 

and fjL 

First we study the Green function for a complex scalar field in the constant 
magnetic field, at finite temperature, T = 1/(3, and chemical potential, /x. 
The chemical potential is defined for the global U(l) symmetry of the complex 
field. In the constant magnetic field, the Green function, G{x,y;ms), for 
the scalar field obeys the Klein-Gordon equation 

{{d^ + ieA^y + ml} G{x,y;ms) = 6\x-y). (1) 

We introduce the temperature and the chemical potential to this equation. 
Following the standard procedure of the imaginary time formalism,^^'^^ the 
thermal Green function is defined in Euclidean space-time by 

[{idi - i/i + eAif + {idj + eAjf + ml] G{x,y;ms) = 5%{x - y), (2) 

where S^{x — y) = 5{xi — yi)5^{x — y). The time direction, X4, is restricted 
between and /?, i.e. X4 G [0,/?]. Thus the thermal theory has no Lorentz 
invariance. It naturally follows that the thermal equilibrium is defined along 
a specific time direction. 

The Klein-Gordon equation is simplified by expanding the time direction, 
Xi, in Fourier series as 

^ 00 

G{^^y) = 3!! e-'^-^'''-'''>Gn{x,y), (3) 

n=— 00 

where the Matsubara frequency Un is given by u;„ = 27m/ (3 for a scalar field. 

Here we consider the constant magnetic field along the z-axis, = 
Sfj^2XiH, for simplicity. For the constant magnetic field the fourth component 
of vanishes, A4 = 0. For a constant Eq.Q reduces to 

{{un - z/i)' - {dj - ieAjf + ml] Gn{x,y) = 6%x - y). (4) 

The induced Green function, G^, has similar form to the three dimensional 
Green function for the scalar field with mass M = a/ {uJn — ^/^)^ + in 
the electromagnetic field. M develops an imaginary part only if both the 
temperature and the chemical potential have non-vanishing value. Because 
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of this imaginary part we must modify the original Fock-Schwinger method 
as is shown below. 

We consider the proper-time Hamiltonian Hnf'^'^ 

3 

Hn = - - leAj)^ + {ujn - ifi^ + ml (5) 

i=i 

Evolution of the system is described by the proper-time, r. Then the induced 
Green function Gn in Eq.(j3]) satisfies 

HnGn{x,y)=6'{x-y). (6) 

To solve this equation we introduce the unitary evolution operators and 
which are defined by 

^^U:{x, y- t) = H^U:{x, y; r), {a = 1, 2), (7) 

with the boundary conditions 

lim Ul{x,y-T) = Urn U^{x,y;T) = 0, (8) 

T — > — OC T — :*00 

lim Ul{x,y;T) = lim U^{x,y;T) = 6'\x - y). (9) 

In the case of vanishing temperature or vanishing chemical potential the in- 
duced Green function, G„, can be described by only one evolution operator, 
U^. However, two types of the evolution operators with different boundary 
conditions are necessary to obtain a finite G„ under non-vanishing temper- 
ature and chemical potential. In the previous works contributions from the 
evolution operator f/^ is not completely considered. The induced Green func- 
tions Gn{x, y) are expressed by U!^ as 

{-i ! dr U^{x,y;T), ( n < ), 
(10) 
J ^n(a;,?/;r), ( n > ). 

After some straightforward calculations, see for example Refs. 17 and 18, we 
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obtain the evolution operators 



eHr 



(47r)3/2 |-7-|3/2 sin {eHr) 



exp {le di- A{i) 



X exp 



{x - y)ieFij[coih. (eFr)]jfc(x - y)^ 



-ir {{un-iiif + ml] , (11) 
where F is the field strength and H is the magnetic field, a"" is defined by 

(12) 



g+37ri/4 
„-37rj/4 



a = 1 ), 
a = 2 ). 



The evolution operators, [/^, are exponentially suppressed at the limit r — ^ 
—oo for > and r ^ cxd for n < 0. It is clearly seen in Eq. (fTT|) that the 
evolution operators obey the boundary conditions (jH)) and ©• Substitut- 
ing Eqs. (fTU|) and (fTTj) in Eq.(j2I) we get the explicit expression of the Green 
function, 



G(x,y]ms 



—I 

T 



n=-l 



-iuJn{x4—y4) 



n=l 



rfrf/^(r) 



-itJn (2:4—2/4) 



dr UI{t) 



le 



■i-Ki/A 



(47r)3/2/5 
dr 



,iaj„ (2:4-4/4) 



n=l 

eH 



^1/2 gjj-^ (^eHr) 



exp / rfe ■ A(0 



X exp 



--{x - y)ieFij[coth {eFT)]jk{x - y)k 



+iT { {ujn + ifi)"^ + ml] 



+ c.c. . 



(13) 



The first term on the final result in Eq. ()13|) comes from the integration of 
the evolution operator, U^. The complex conjugate of this, (c.c), appears 
from which is introduced to deal with the non- vanishing temperature and 
chemical potential. 
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3 Spinor Two-point Function at Finite T 

and 11 

Next we consider the Green function for a fermion field in an external con- 
stant magnetic field at finite temperature and chemical potential. The Green 
function, S{x,y;'mf), is defined by the Dirac equation; 

{i^ + e^.- ifi-f4- TUf) S{x,y;mf) =6^{x-y), (14) 

where m/ is the mass of the fermion field. 

To calculate the analytical form of S{x,y]mf) it is more convenient to 
introduce the bi-spinor function Gf{x,y), 

S{x, y;mf) = {i^ + — il^li + ^f) Gf{x, y). (15) 

The explicit form of S{x, y; rrif) is determined by solving the following equa- 
tion for Gf{x, y), 

+ ^ei/(7i72 - 727i) - 2fxd, + /x^ - mjj G/(x, y) = 4(x - y), (16) 

where Dj = dj —ieAj. Here we choose the direction of the constant magnetic 
field along the z-axis, = 5^2XiH. For a constant magnetic field the 
function, Gf{x,y), is expanded in Fourier series 

oo 

Gf{x,y) = - e-'^-^^^-y^^Gi{x,y), (17) 

n=— oo 

and Eq. (fT^ reads 

jS^i + ^ei/(7i72 - 727i) - (^n - - mjj Gi{x,y) = S%x - y). 

(18) 

As in the scalar case, the induced function, Gl{x,y), is calculated by intro- 
ducing the proper-time Hamiltonian, 

3 

Hi = J2d' + ^ei^(7i72 - 727i) - (^« - ^/i)' - mj. (19) 
i=i 
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According to the similar way in the scalar field we can find the bi-spinor 
function, Gf{x,y). It is described by two types of the unitary evolution 
operators and which satisfy the boundary conditions (jHl) and 
Therefore the explicit expression of Gf{x,y) is obtained by 



Gf{x,y) 



-oo 



n=0 



-oo 



-E 



-iuJn{x4-y4) 



n=-l 



POO 

J+o 



■ (20) 



The Matsubara frequency for the fermion field is given by u;„ = {2n-\- 1)ti / (3. 
The evolution operators U^{t) are found to be 



Kir) 



eHr 



(47r)3/2|r|3/2 sin(effr) 



exp j di ■ A{i) 



X exp 



--(x - y)ieFij[coth {eFT)]jk{x - y)k 



-IT < -eFjkCTjk - {uJn - if^y - mj 



(21) 



where ajk = -[7j,7fc], F is the field strength and 6" is 



6" 



3-37ri/4 
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3+37ri/4 



a = 1 ), 
a = 2). 



(22) 



Inserting Eq. ()2ip into Eq. lj^UI) . one easily derive the two-point Green function, 

S{x,y;mf) = 



eH 



le 



-3m /A \ 



(47r)3/2/? 



^ n=0 



-iuj„(x4— 1/4) 



POO 

X dr 
Jo 

X exp 



rV2sin(e/7rT'"Pi"y/^-^^^) 



(x - ?/)ieFij[coth {eFT)]jk{x - y)^ 



+ (c.c). 



(23) 



The complex conjugate term, (c.c), in Eq. (j23j) comes from which is intro- 
duced to deal with effects of both the temperature and the chemical potential. 

4 Behavior of the Vacuum Self-energy 

At the one loop level the vacuum self energy for a free field is given by the 
trace of the Green function. Here we calculate it for a free scalar and a free 
fermion fields at finite H, T and fi. 

For a free scalar field with mass the trace of the Green function ()13p 
becomes 



where V is the 3-dimensional volume. If it is larger than m^, naive perturba- 
tion loses validity to evaluate the radiative correction in a scalar theory with 
interactions. To get rid of this difficulty we must use a resumed propagator 
known as the ring diagram resummation in the thermal field theory. 

Carrying out the trace of the Green function (PHj) with respect to space- 
time and spinor legs, we find the vacuum self-energy for a free fermion field, 



For the vanishing temperature and/or chemical potential the contribution 
from the evolution operator is coincide with the one from U^^ and contour 
dependence of the proper-time integral disappears. Therefore our result must 
be coincide with the one found in the previous work. For /i — Eq. ()25p 
exactly agrees with the one obtained in Ref. 6. At the limit, T ^ (/? — > oo), 
Eq.((2ni) reproduce the results obtained in Refs. 8 and 9. 

Performing the proper-time integral and the summation in Eq. ()24|l and 
Eq. ()25p numerically, we evaluate the vacuum self-energy for the free scalar 
and the free fermion fields. We choose the contour of proper-time integration 
in the first term of the right hand side of Eq. ()24|) slightly above the real axis. 
As is known, this contour gives physical results at the limit /i ^ and 
T 0.^ It is natural to take the same contour at finite fi and T. 




X exp [zr {(u;„ + ifif + m^}] + (c.c), (24) 




X exp [-ir {(cj^ - ifi^ + mj]] + {c.c). (25) 



8 




poles 

(a) (b) 
Figure 1: Contour on the complex r plane. 



For uf — fi"^ + ml > the integrand is exponentially suppressed at the 
infinity above the real axis. We close the contour as is shown in Fig. |3] (a) and 
perform the proper-time integral along the path Ci. The complex conjugate 
of the result gives the second term in the right hand side in Eq. ()24|) . Thus 
the TiG{x, x) is found to be 



(3V 

where /s(t, n) is 



1 1 _ /"OO 

X) = J] j^^^^ drUr, n), (26) 



/.(,.„) = £^ "=P-"^"l -^»".'.-^'^-;). (27) 
a/t smh(eiir) 

Here we introduce the proper-time cut-off A to regularize the theory. 

In the case uf — fi"^ + mj. < we must consider two kinds of paths in 
Fig. m (a) and (b). For a positive tu^ — /i^ + the integrand drops at the 
infinity above the real axis. If — fi"^ + is negative, the integrand is 
suppressed at the infinity below the real axis. We calculate the proper-time 
integral along the paths Ci, C2, C3 and add the contribution from the poles 
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on the real axis. After some calculations TrG{x, x) is obtained by 



^TtG{x, x) 



^ oo 
V2^ 



47r3/2/5 



n>[N] 
[N] 



dTfs{T,n) 



1/A2 



3/2/3 



47r3/2/5 ^ 

71=1 '- 



dTgs{T,n) 



1/A2 



,(28) 



where = (3^y fi'^ — m2/(27r), [A^] is the Gauss notation and 



/iso(^) 



eH sin(2a;n/xr) „^ 
r sinh(e-ffr) 



_^/2 A sin(eii'e^7A2) 

21 „ie / A 2i 



X exp[z{(cj„ + z/i)^ + }e'7A^] + (c.c.) 



(29) 



(30) 



oo /■ 



1/2 



-2iTlu>ntJ./ (eH) 



X COS 



/ 2 2 2\ ^ 

— (/i -^„-mJ + - 



(31) 



hso{n) is the contribution from the path C3 around the pole at r = and 
hsj{n) is sum of residues at eHr = ttI. 

In the case of the fermion field the above contour gives the physical re- 
sult at yU and T — »• for the second term "(c.c.)" of right hand side of 
Eq.(j2Hl). According to the similar analysis with the scalar field, we calcu- 
late the TyS{x,x). For a positive lOq — fi'^ + ml the proper-time integral is 
performed along the path Ci, 



1 

where f{T,n) is 



TyS{x, x) 



eH 



t3/2/3 



7r3/2/5 



n=0 



dTf{T,n) 



(32) 



1/A2 



/(r,n) = ^coth(ei/r) cos(2cj„/ir)e-"(""-^ 

0" 



(33) 
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For a negative cJq — /i^ + we evaluate the proper-time integral along the 
path Ci, C2, C3 and consider the influence from pole on the real axis. 



TtS{x, x) 



TTli 



ruf 



V / dTf{T,n) 

n>[N] -^l/A^ 
[N] 



n=0 



ho{n) + hj{n) + / dTg{T,n) 

1/A2 



(34) 



where N 



fi'^ — mj/(27r) — 1/2, g{T,n), ho{n) and /ij(T2) are given by 



hoin) 
hj{n) 



cothieHr) sin(2cu„/ir)e-"(^' 



37rj/4 rTr/2 tt 

/ rf^£fie^^/2^ot(e//e^7A2) 
X exp[z{(u;„ + i^if + mJje'VA^] + (c.c.) 



(35) 



(36) 



X sm 



tt/ 



(37) 



/io corresponds to the contribution from the path C3 and hj is sum of residues 
at eiir = nl. 

Performing the integration over r and the summation numerically, we 
obtain the behaviors of the trace of the Green functions, i.e. the vacuum self 
energy at the one loop level. In Fig. |21 we show behaviors of the vacuum 
self-energy as a function of the external magnetic field H with T and /i fixed. 
For a case of neutron star many interest has been payed to the QCD phase 
structure. Thus we suppose the proper-time cut-off is more than QCD scale, 
^QCD < A ~ 2GeV. The scalar and fermion mass is taken to be the pion 
mass scale, ~ ~ O.lGeV. The present upper limit for the magnetic 
field in the neutron star is of the order, eH ~ O(O.OlGeV^).^^ An oscillating 
mode is observed for both the scalar and fermion field. The amplitude of the 
oscillation becomes larger as H increases and/or T decreases. The oscillation 
disappears for higher temperature. For the neutron star the upper limit of 
the magnetic field, eH, is of the order O(0.01GeV~^). It seems to be difficult 
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Figure 2: Behaviors of the vacuum self energy as a function of H with T 
and /i fixed. We set A = 2GeV, = mj = O.lGeV, fi = IGeV and 
T = 0.006GeV, T = O.OOSGeV, T = O.OlGeV. 




Figure 3: Behaviors of the vacuum self energy as a function of /i with H 
and T fixed. We set A = 2GeV, m, = m/ = O.lGeV, T = 0.006GeV and 
eH = O.OOOlGeV^, eH = 0.02GeV2, eH = 0.04GeV2, eH = O.OeGeV^. 
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to see the magnetic oscillation appeared in Fig. |21 unless the neutron star is 
extremely cold. It agrees with the result obtained in Ref. 11. 

Behaviors of the vacuum self-energy is illustrated as a function of fi with 
T and H fixed in Fig. El The trace of the two-point function goes down as 
H increases for both the scalar and fermion field. We can see the oscillating 
mode for both cases. For the fermion field the mode is the origin of the van 
Alphen-de Haas magnetic oscillations as is shown in Ref. 11. Such a effect is 
not found in the scalar field theory, since the scalar field has no sharp Fermi 
surface. 

As in known, the trace of the scalar two-point function contains a term 
which is proportional to T"^}^ Indeed, Eq. (l27|) reduces to the well-known 
result, 

^-2- + 0(7). (38) 

at the limit fi, H and — > 0. To obtain it we drop the surface term of 
the proper-time integral and use a formula C(~l) = —1/12. behavior 
is not observed in Fig. |21 because the surface term is proportional to AT. 
The temperature considered here is too small compared with the cut-off 
scale A. T-dependence of the surface term is canceled out if we take the 
T-dependent cut-off, A oc 1/T. Similar property is found in the momentum 
cut-off regularization for H = 0. 



5 Conclusion 

We have investigated the scalar and fermion field theories at finite tempera- 
ture, chemical potential and constant magnetic field. The explicit expressions 
of the two-point Green functions are found by using the proper-time formal- 
ism. If both the temperature and the chemical potential exist, we must 
modify the Fock-Schwinger proper-time method by introducing two types of 
the evolution operators with different boundary conditions. 

The proper-time integrations remain in our final expressions of the two- 
point Green functions G{x, y; rris) and S{x, y\mf). The remained integrand is 
exponentially suppressed at the limit r ^ oo. There are poles at cHt = n% 
for any integer n. Because of these poles the naive Wick rotation has no 
validity. ^ We carefully take the contour of proper-time integration slightly 

^ After the naive Wick rotation t it, we can perform the summations in Ea. H24l) and 
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above the real axis to avoid these poles in the complex r plane. In the case 
with large chemical potential we must consider two kinds of contour. One of 
the contour contains the poles at eHr = nir. On the other hand, the naive 
analytic continuation from the low chemical potential takes the contour below 
the real axis for Un — fJ-"^ + rn? < 0. Thus it drops the contribution from the 
poles on the real axis. It gives only approximate results. 

We apply our formalism to the vacuum self-energy at the one-loop level. 
At the limit T and/or /i — >■ contributions from the poles at eHr = n% 
disappear and our results coincide with the previous one. We numerically 
perform the proper-time integral and the Matsubara mode summation. In 
a strong magnetic field oscillating mode is observed for both the scalar and 
the fermion fields. Our results qualitatively agree with the previous analysis 
obtained by the naive Wick rotation^^ for fi < rrif. Performing the naive 
Wick rotation, Eq.(j2S)) reads 

— TrS = -— ^ / dr^ coth(ei/r)e-(^'-'"?)^2(2/ir//3,4vr^r//32) . 

(39) 

In Fig. I^we show the behaviors of the fermion vacuum self-energy p9|l . T 
and /i are fixed on the same value with the Fig. IHl We cut the proper-time 
integration at lO^GeV"^ because it does not converge for ^ > rrif. As is 
shown in Fig. HI the vacuum self-energy coincides with the results in Fig. El 
for fi < TUf = O.lGeV. But it is completely different for /i > O.lGeV. In Fig. 
4 the vacuum self-energy is not analytic at some points. 

There are some interesting applications of our work. Using the two-point 
functions obtained here we calculate the radiative correction of some QED 
process, radiative decay of axion, neutrino, and so on. As an example, we 
calculated the trace of the two-point functions. It is necessary to determine 
the phase structure of the theory at finite T, /i and H."^^ The influence of 
the external magnetic field will modify the decay rate and phase structure 
at finite T and fi. It may affect the evolution of our universe and/or neutron 
stars. 

The present work is restricted to the analysis of the influence of the 
temperature, chemical potential and a constant magnetic field. There are 
some interesting objects in an external electromagnetic field. However it is 

Ea. (P3l) and the vacuum self-energy is described by the eUiptic theta function, However, 
it is not vahd in our formalism. 
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Figure 4: Behaviors of the vacuum self energy as a function of /i with H 
and T fixed. We set A = 2GeV, = m/ = O.lGeV, T = 0.006GeV and 
eH = O.OOOlGeV^, eH = 0.02GeV2, eH = 0.04GeV2, eH = O.OeGeV^. A 
vertical line of right side consists of all the value of eH. 



not clear to extend our procedure to the state in an external electric field. It 
is also interesting to introduce the curvature effects in our analysis. 
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